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Abstract 

! In this article we evaluate mass corrections up to 0{(m? /g 2 ) 6 ) to the three-loop 

polarization function induced by an axial- vector current. Special emphasis is put 
on the evaluation of the singlet diagram which is absent in the vector case. As a 
physical application 0{a 2 s ) corrections to the production of top quarks at future 
e + e~ colliders is considered. It is demonstrated that for center of mass energies 
y^s ;> 500 GeV the inclusion of the first seven terms into the cross section leads to 
a reliable description. 

X ' PACS numbers: 12.38.-t, 12.38.Bx, 13.85.Lg, 14.65.Ha. 

a: 

In the total cross section cx(e e — > hadrons) corrections arising from the finite mass, 
m, of the produced quarks may often be neglected. Concerning precision measurements 
around the Z resonance first order mass corrections, known up to 0(a^) [p], ||, are usually 
adequate. However, having in mind top quark production at future colliders like the NLC 
with a center of mass energy of \/s = 500 GeV higher order terms in m? / s may become 
important. The velocity of the produced particles is then v ~ 0.7 which means that on 
one side threshold effects are not important and on the other side we are not in the region 
of very high energies. 

In Ref. Q] the contribution of the photon to the production of top quarks was con- 
sidered. In this article also the exchange of the Z boson is included. Hence, in a first 
step results for the axial-vector polarization function up to 0(a 2 s ) are presented. The 
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imaginary part in combination with the recently evaluated rate for the vector case 
directly leads to the cross section cr(e + e~ — > tt + X) mediated by a virtual Z boson. The 
0{a s ) corrections to this process were considered in ||. 

To be more precise let us define the axial-vector current correlator as: 



(-<? V + Qtfv) n V) + w* n 2(? 2 



i dxe^(0\Tf(x)jZ(0)\0) 
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with = ip^^y^. In the following we will only present results for Il a (g 2 )[|. It is convenient 
to write 
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with the SU(3) colour factors Cp = 4/3, Ca = 3 and T = 1/2. Yi^ ,a is the abelian 
contribution already present in QED and II^ a originates from the non-abelian structure 
specific for QCD. The polarization functions containing a second massless or massive 
quark loop are denoted by n| 2 ^ a and IT}/ > ' a , respectively. IP represents the so-called non- 
singlet part. However, for external axial-vector currents already at 0(a 2 ) there exists also 
a singlet or double-triangle contribution: 
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As IT| depends on the properties of both members of the fermion doublet we will from now 
on specify to the top-bottom case. The generalization to other quark flavours is obvious. 
For this contribution it is convenient to replace the current j* in Eq. ([[]) by tj^t — bj^b 
because in this combination the axial anomaly cancels. In Fig. |I| the relevant diagrams 
are depicted. 

Similar relations as in Eqs. (|2]) and (|^) also hold for R a (s) and Rg(s), respectively, 
defined through 



= 127rlmU a {s) (q 2 = s + ie), (4) 
so that the cross section for the inclusive production of top quarks may be written as 
a(e + e- -^tt + X) 



Rt(s) 



a 



pt 



(v 2 e + a 2 e )a 2 t 



s-M 2 Z/ 



R a (s) + R a s (s) 



71 



— CpTR 



Sb \ 



QiQf + 2Q e v e Q t v t s _ M2 + (v 2 e + 4)v 



2\„.2 
t 



R v (s) 



(5) 



1 The longitudinal part, H1(q 2 ), of the non-singlet contribution, e.g., is via the axial Ward identity 
directly connected to the pseudo-scalar polarization function W(q ), for which the high energy expansion 
was considered in 101. 
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Figure 1: Diagrams contributing to Ug. In the triangle loops either a top or bottom quark 
may be present. 



with a pt = Ana 2 /3s, v f = (// - 2Q f s$)/(2s 9 c 9 ), a f = ll/(2s e c e ), Q e = -l,Q t = 2/3, If = 
— 1/2 and i| = 1/2. Furthermore we have c 2 , = 1 — s 2 , with sg being the sine of the weak 
mixing angle. R v (s) is given in and both R a (s) and will be presented below. 

R^Sb is the contribution from cuts of the singlet diagram that do not involve top quarks. 
Non-singlet contributions with the photon or Z boson coupling to a light quark flavour 
and the top quarks produced via gluon splitting |7| will be neglected as their numerical 
values are tiny [||. 

The computation of IP naturally splits into two parts: Firstly into the non-singlet 
contribution where the anticommuting definition of 75 may be used. Here the calculation 
of the diagrams is in close analogy to the vector case. Hence we refer for details to [|J. 

The second part, the singlet contribution IT|, is connected with the axial anomaly and 
is not present in II". Let us briefly describe our treatment of these diagrams. Actually 
three graphs have to be considered, namely the cases when two top quarks, one top and 
one bottom quark or two bottom quarks are running in the triangle loops. One may argue 
that the last combination only contributes to the cross section into bottom quarks which 
is not the process under consideration. However, only the proper combination of all three 
parts guarantees the cancellation of the anomaly. From the final result the cuts arising 
from bottom quarks have to be subtracted, of course. 

For the evaluation of n^' a naive 75 fails to work. We follow the treatment introduced 
in || and formalized in || and replace both axial-vector vertices according to [K| 

7„7s - 7^A PCT 7 [AH , (6) 

where ^ xpa ^ is the antisymmetric combination of three 7 matrices which can be written 
as ^,I A H — ^ x ^p^°~ _ 'j a 'jp'j x y2. In a first step the e-tensors are put aside and the new 
object with six external indices, n^y CT ,j, defined through 
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is treated until the momentum integration and renormalization is done and a finite quan- 



tity is available |11J. Then the contraction with the e-tensors is performed. It is possible to 
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show that the contribution from the singlet diagrams may be computed from the relation 
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which means that we can treat the scalar quantity g CT g^Il|j^L in complete analogy to 
the non-singlet diagrams. We should mention that a finite renormalization of the singlet 



axial-vector current |13[ has not to be performed in the order considered in this paper. 

Using the large momentum procedure the first seven terms in the m 2 /g 2 -expansion 
of Il a (g 2 ) have been evaluated. We refrain from listing the results separated into the 
contributions from the different colour factors and present the results for the proper sum 
keeping only ni, the number of light (massless) quarks, as arbitrary parameter (l qm = 
ln(-g 2 /m 2 ),/ g ^ln(-g 2 /^ 2 )): 
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where m t is the MS top mass and ( is Riemann's zeta-function with the values (2 — k 2 /6, 
C 3 « 1.20206, (4 = tt 4 /90 and (5 ~ 1.03693. B 4 « -1.76280 is a constant typical for 
massive three-loop integrals ||14|| . The expansion of the two-loop quantity, n^' can be 



compared with the exact result |15[ . At order a s the constant and quadratic terms are in 
agreement with [[U], |17|. Note that in the non-singlet contribution m t could be replaced 
by any other quark mass. The singlet part, however, gets modified if both quarks have 
to be considered as massive and even vanishes for a degenerate quark doublet. This is 
also the reason for the absence of the first two terms in the expansion for m 2 /q 2 — > 0: 
For m t = the top and bottom quark are trivially degenerate and the contribution to 
the first order power corrections arise from a simple expansion of the diagrams for small 
masses. According to the structure of the 7 matrices from each triangle at least a factor 
mf has to come. This means that the m 2 / q 2 corrections from the diagrams with two top 
triangles cancel against the one with a top and a bottom triangle which has an overall 
factor of two. 

Taking the imaginary part of Eqs. (|9|-|r2|) and transforming the result into the on-shell 
scheme concerning the top mass leads to (L ms = ln(M 4 2 /s)): 
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where /z 2 = s is chosen. Note, that the quartic corrections of IIg have no imaginary 
parts so that R^' a actually starts at order (M t 2 /s) 3 . An important check of our result is 
provided by the successful comparison of the terms proportional to ni with the expansion 
of the exact analytical expression [R3|. The quartic terms for the proper sum are 
also available in the literature p0| and complete agreement was found. 

For completeness we list the results from the double-triangle diagrams containing cuts 
from the b quark only pi] : 
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This contribution has to be subtracted from Rg' a . Note that the cut arising from two 



gluons is zero according to the Landau- Yang-Theorem p^| . In Fig. |^ the terms for the five 
different contributions are plotted against x = 2M t /y/s including successively higher or- 
ders in Mils. For i$ ),a and R$f a comparison with a recently evaluated semi-analytical 
result (narrow dots) [O is possible and agreement up to x ~ 0.8 is found. The light 
fermion contribution, i?[ 2 , may be compared with exact results JTP| (narrow dots) and 
also shows agreement up to x ~ 0.8. Concerning R^' a and R^ ,a the situation is less 
satisfactory. It seems that there is reasonable convergence up to x ~ 0.7 which is also 
motivated by the behaviour of the vector case where analytical results for x > 0.5 are 
available (see ||). However, for x > 0.7 the behaviour of the curve including all known 
power correction terms (solid line) indicates that close to x — 1 the convergence fails to 
work. The reason presumably is connected to the four particle cut starting at x — 0.5. 
Although R^' a and R$a also exhibit a four particle cut it seems to be somehow less 
dominant for these contributions. In Fig. ^| also the contribution R^' a is shown. Here, 
already the curve including power corrections up to order (M 2 /s) 2 is practically indis- 
tinguishable from the exact result. Note that only the difference of the two plots in the 
bottom line of Fig. Renters Rt{s). 

Recalling that the first seven terms for R v (s) approximate the exact result up to 
x « 0.7 |4]] we are now prepared to present predictions for R t valid up to 0(a 2 s ) for 
y/s J> 500 GeV which corresponds to x <> 0.7. Therefore we insert the isospin and charge 



quantum numbers into Eq. @ and choose n\ = 5, 
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Mz = 91.187 GeV and M t = 175 GeV. Then the expansion of Rt(s) looks as follows: 
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In Tab. |T| the coefficients R^ are listed for different values of the center of mass energy 
\fs. One observes that for y/s = 500 GeV, which is a proposed option for the NLC, the 
O(ol) QCD corrections amount to ~ 2%. For higher values of the center of mass energy 
these terms get less important. In Fig. ||] the normalized cross section R t is plotted against 
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Figure 2: R\ 2) ' a , i = A, NA, I, F , S , Sb as functions of x = 2M t / ' y/s at fi 2 = Mf. 
Successively higher order terms in (M t 2 /s) n : Dotted: n = 0; dashed: n — 1, ... ,5; solid: 
n — 6. Narrow dots: exact result {R^' a , Rsb' a ) or semi-analytical results (R$' a , R^na)- 



y/s. The contributions from the vector and axial- vector part are also displayed separately. 
Rt(s) is clearly dominated by the vector contribution which is mainly due to the fact that 
in Eq. (5) the couplings to R v are larger by roughly a factor of four as compared to R a . 
Another reason is that the Born cross section R^' v is always larger than R^' a . This is 
not true for the 0(a s ) and 0(a 2 s ) terms. Here the axial-vector contribution exceeds the 
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R? ] 


Rt(s) 


500 


0.70 


0.095 


1.419 


6.021 


29.902 


1.629 


1000 


0.35 


0.088 


1.732 


2.842 


6.016 


1.816 


1500 


0.23 


0.085 


1.771 


2.291 


3.709 


1.836 


2000 


0.18 


0.083 


1.784 


2.091 


3.044 


1.841 



Table 1: Numerical values for the contributions of 0{a\), R t , to the normalized cross 
section R t . The values for af\s) are based on a^(M|) = 0.118. The scale fj 2 = s has 
been adopted. Also the values of x = 2M t / \fs are shown. 

2 
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1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 

° 600 800 1000 1200 1400 1600 1800 2000 

Vs~[GeV] 

Figure 3: Normalized cross section R t as a function of the center of mass 
energy y/s, together with the pure vector and axial- vector contributions: R t = 
^vector + ^axial-vector _ Dottec i; Born approximation; dashed: 0(a 8 ), solid: 
0(a 2 s ). The scale fi 2 = s has been adopted. 



vector part for sufficiently large values of \fs and approaches it from above as \fs goes to 
infinity, where both R v and R a are identical. In Fig. |] this is demonstrated at order a 2 . 
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Figure 4: Two-loop vector and axial- vector contributions R^' v , R^' a + 
CpTRg' a . The scale fj 2 = s has been adopted and M t = 175 GeV has been 
chosen. 

For this reason at energies above roughly ^Js = 600 GeV the tree-loop vector and axial- 
vector contributions to R t are comparable. At lower values of the energy the vector part 
is still larger than the axial- vector part as a consequence of the more singular threshold 
behaviour of R v (Fig. |4|). For the sake of completeness we note that the contribution from 
the singlet diagram which is absent in the vector case is smaller by at least a factor 100 
as compared to the non-singlet case. 

To conclude, the large momentum procedure has been applied to the axial-vector 
polarization function and terms up to order (M 2 /q 2 ) 6 have been determined. The imag- 
inary part in combination with the result recently obtained for the vector case was used 
to predict the production of top quarks at future e + e~ colliders up to 0(a 2 ). 
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